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Dark matter particles annihilating in the Galactic halo can provide a flux of positrons potentially 
observable in upcoming experiments, such as PAMELA and AMS-02. We discuss the spectral 
features which may be associated with dark matter annihilation in the positron spectrum and assess 
the prospects for observing such features in future experiments. Although we focus on some specific 
dark matter candidates, neutralinos and Kaluza-Klein states, we carry out our study in a model 
independent fashion. We also revisit the positron spectrum observed by HEAT. 
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I. INTRODUCTION 

The existence of dark matter has been confirmed by a 
wide array of experimental tests including observations 
of galactic clusters and large scale structure [1], super- 
novae [2] and the Cosmic Microwave Background (CMB) 
anisotropies [3,4]. WMAP has measured the density of 



0.113±g:8i| [4] at the 



cold dark matter to be ^cdm^ 2 
2a confidence level. 

For a variety of reasons, it is often thought that dark 
matter may be made up of weakly interacting, TeV-scale 
particles [5]. If this is the case, there are a variety of 
experimental approaches which may be capable of ob- 
serving them directly or indirectly. The most straight 
forward approach is to observe dark matter particles di- 
rectly as they scatter off of a detector [6] . Such particles 
could also be produced and observed in colliders such as 
the Tevatron or the Large Hadron Collider (LHC) [7]. Al- 
ternatively, dark matter particles which annihilate in our 
Galaxy's halo or in gravitational wells such as the Galac- 
tic Center or the Sun may produce potentially observ- 
able fluxes of gamma-rays [8] , neutrinos [9] , anti-protons 
[10,11], anti-deuterons [11,12] or positrons [11,13-17]. 

The HEAT (High-Energy Antimatter Telescope) ex- 
periment, in three flights taking place in 1994, 1995 and 
2000, observed a flux of cosmic positrons in excess of 
the predicted rate, peaking around 10 GeV [18-20]. Al- 
though the source of these positrons is not known, it has 
been suggested in numerous publications that this sig- 
nal could be the product of dark matter annihilations 
[13-17]. 

Although HEAT provided the most accurate mea- 
surement of the cosmic positron spectrum at the time, 
these flights had their limitations. Due to the rapidly 
falling flux, HEAT was able to measure the spectrum of 
positrons only up to approximately 30 GeV, and with 
rather large error bars down to 15 GeV, or so. Fortu- 
nately, other experiments, such as PAMELA and AMS- 
02, are planned to improve considerably on these mea- 
surements. 

PAMELA (Payload for Antimatter-Matter Explo- 



ration and Light-nuclei Astrophysics) is a satellite borne 
experiment designed to study the matter-antimatter 
asymmetry of the universe with very precise cosmic ray 
measurements. PAMELA'S primary objectives include 
the measurement of the cosmic positron spectrum up to 
270 GeV [21], well beyond the range studied by HEAT. 
This improvement is made possible by PAMELA'S large 
acceptance (20.5 cm 2 sr [22]) and long exposure time (3 
years). PAMELA is scheduled for launch in 2005. 

AMS-02 (Alpha Magnetic Spectrometer) is an experi- 
ment designed to be deployed on the International Space 
Station (ISS) for a three year mission sometime around 
the end of the decade, perhaps 2008. AMS-02 's accep- 
tance of about 450 cm 2 sr is considerably larger than 
PAMELA'S. This, along with superior energy resolu- 
tion and electron and anti-proton rejection, makes AMS- 
02 the premiere experiment for measuring the cosmic 
positron spectrum. 

The precision measurements of the cosmic positron 
spectrum to be provided by PAMELA and AMS-02 will 
be of critical importance to efforts to identify signatures 
of dark matter in this channel. Although the overall 
flux of positrons is important, without detailed spectral 
information it will be difficult to distinguish the prod- 
ucts of dark matter annihilation from other possible as- 
trophysical positron sources. Possible sources of cosmic 
positrons include their production via hadronic cosmic 
ray interactions with giant molecular clouds or electron- 
positron pair creation via electromagnetic interactions in 
nearby pulsars. Alternatively, radioactive nuclei ejected 
in supernova blasts could potentially produce high energy 
positrons [20]. 

In this article, we discuss the spectral characteristics 
of cosmic positrons associated with the annihilation of 
dark matter particles in the Galactic halo and assess the 
prospects for their detection by PAMELA and AMS-02. 
The remainder of this article is organized as follows. In 
section II, we discuss the spectrum of positrons produced 
in particle dark matter annihilations and discuss the ef- 
fects of propagation on this spectrum. In section III, we 
show the positron spectra which result for several dark 
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matter annihilation channels. In section IV, we discuss 
possible variations to these results from our choice of dif- 
fusion parameters and dark matter halo profile. In sec- 
tions V and VI, we describe the positron production of 
two specific dark matter candidates: neutralinos in mod- 
els of supersymmetry and stable Kaluza-Klein states in 
models of Universal Extra Dimensions (UED). In section 
VII, we discuss our results in the context of the measure- 
ments made by the HEAT experiment. In sections VIII 
and IX, we discuss the prospects for the detection of dark 
matter annihilation in the future experiments PAMELA 
and AMS-02. In section X, we compare this technique 
to other dark matter detection methods. Our concluding 
remarks are contained in section XI. 



II. POSITRON PRODUCTION AND 
PROPAGATION 

Positrons can be produced in the annihilations of dark 
matter particles through many channels. Annihilations 
which yield gauge bosons, for example, yield positrons 
as these gauge bosons decay, peaking at an energy near 
<~ mj/2. In addition to the decays Z — > e + e~ and 
W + — > e + v, positrons can be produced in gauge bo- 
son decays to muons, Z — ► yU + /i~ — > e + vve~vv or 
W + — > n + v — > e + vvv. Gauge bosons which decay to 
tau pairs can produce positrons directly in their decays, 
through their decays to muons, or hadronically through 
the decay of charged pions. Of course, gauge bosons may 
also decay to hadrons which also produce positrons via 
charged pions. 

Dark matter annihilations to bottom quark pairs are 
common, especially in the case of bino-like, neutralino 
dark matter. These b quarks decay, producing charged 
leptons without the feature seen in gauge boson decays 
near <~ Well below this energy, the positron spec- 

trum from b decays is similar to that for gauge boson 
decays. 

Dark matter annihilations to top quarks always pro- 
duce a pair of bottom quarks and W ± bosons, which can 
each generate positrons as described above. As would be 
expected, annihilations directly to charged leptons are ca- 
pable of producing a harder spectrum of positrons than 
the other modes one could consider. In figure 1, we com- 
pare the positron spectrum from particle annihilations 
via several channels. Annihilations to gauge bosons or 
heavy quarks each produce a much softer spectrum than 
the t + t~ or n + channels. Although it is not shown, 
annihilations directly to e+e~ would be represented sim- 
ply by a delta function at an energy equal to the WIMP 
mass. For the results shown in figure 1, and throughout 
this paper, the spectrum from each annihilation channel, 
including cascading, is calculated using P YTHIA [23] , as 
it is implemented in the DarkSusy package [24] . 
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FIG. 1. The positron spectrum from dark matter annihi- 
lations prior to propagation for several annihilation channels. 
Dotted and dashed lines represent the positron spectrum, per 
annihilation, for annihilations to bb and gauge boson pairs, 
respectively. The dot-dashed and solid lines represent an- 
nihilations to t + t~ and n + pT , respectively, which produce 
considerably harder spectra. The spectrum for annihilations 
to e + e~ is not shown, but is simply a delta function at the en- 
ergy equal to the WIMP mass. For all cases shown, a WIMP 
mass of 300 GeV was used. This figure originally appeared in 
Ref. [17]. 

The spectra shown in figure 1 are not those observed at 
Earth, however. As positrons travel through the Galac- 
tic halo, they move under the influence of the tangled 
interstellar magnetic fields and lose energy via inverse 
Compton and synchrotron processes. We can take these 
processes into account by solving the diffusion-loss equa- 
tion: 



d dn„ 



dt dE, 



V 



+ 



dE, 



e+ 



- dn e + 
+ > X)V dET + 
dn e + 



K(E, 
b(E e +,x) 



dE„ 



+ Q(E e+ ,x), (1) 



where dn e + / dE e + is the number density of positrons per 
unit energy, K(E e + , x) is the diffusion constant, b(E e + , x) 
is the rate of energy loss and Q(E e + , x) is the source term. 
We parameterize the diffusion constant [25] by 



K(E e 



3.3 x 10 28 



3 0.47 



0.47 



cm 2 s 



and the energy loss rate by 

b(E e+ ) = 10~ 16 E 2 e+ s- 1 . 



(2) 



(3) 



b(E e +) is the result of inverse Compton scattering on 
both starlight and the cosmic microwave background [26] . 
The diffusion parameters are constrained from analyz- 
ing stable nuclei in cosmic rays (primarily by fitting the 
boron to carbon ratio) [27]. 

In equations 2 and 3, we have dropped the dependence 
on location, treating these as constant within the diffu- 
sion zone. For the boundary conditions of the diffusion 
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zone, we consider a slab of thickness 2L, where L is fit 
to observations to be 4 kpc [25,27]. Beyond the bound- 
aries of our diffusion zone, wc drop the positron density 
to zero (free escape boundary conditions). 

Although the technique used here does not explicitly 
include the effects of re-acceleration, these effects are in- 
troduced though the energy dependence of the diffusion 
constant. The effect of re- acceleration is also not impor- 
tant at energies above a few GeV, where we are primarily 
interested. 

In addition to the propagation effects described above, 
as positrons approach the solar system, their interaction 
with the solar wind and magnctosphere can become im- 
portant. These effects, called solar modulation, can be 
roughly parameterized using the technique of Ref. [28]. 
Alternatively, if one assumes that the effects of solar mod- 
ulation are charge sign independent, their impact can be 
removed by considering the ratio of positrons to positrons 
plus electrons at a given energy rather than the positron 
flux alone. This quantity, called the positron fraction, is 
often used in lieu of the positron flux to minimize the un- 
certainties associated with modelling the impact of solar 
modulation. 



III. THE COSMIC POSITRON SPECTRUM 
FROM DARK MATTER ANNIHILATIONS 

The positron spectra produced in dark matter annihi- 
lations, after taking into account the effects of propaga- 
tion, are shown for various channels in figures 2-5. Many 
of the features described in the previous section can be 
identified here. In particular, the rapid decline in the 
spectrum at around mx/2 can be seen in figure 3, cor- 
responding to the energy at which direct production of 
positrons from gauge boson decays no longer contributes. 
In figure 4, an even more dramatic effect is seen at an 
energy equal to the WIMP mass, corresponding to the 
threshold for the direct production of e + e~. 

In these figures, the effects of solar modulation are not 
included. As discussed in the previous section, we have 
chosen to avoid the uncertainties associated with this ef- 
fect by studying the ratio of positrons to positrons plus 
electrons rather than the positron spectrum alone. In 
figures 6-9, the positron fraction is shown for the same 
annihilation channels. To calculate this quantity, the flux 
of primary and secondary electrons as well as secondary 
positrons must be known. We have used the primary and 
secondary fluxes found in Ref. [29]. Also, in figures 6-9, 
the positron data from the HEAT experiments is shown 
for comparison. 

The magnitude, or normalization, of the positron 
fluxes shown in figures 2-9 are set by the annihilation 
rate of dark matter particles. In addition to scaling with 
the WIMP's annihilation cross section and the square of 



the local dark matter density, the annihilation rate de- 
pends on the degree of inhomogeneity in the local dark 
matter distribution. This effect is described by a param- 
eter called the boost factor, BF, which is defined by 



BF = 



J P 2 dV 



(4) 



where the integral is performed roughly over the volume 
which contributes substantially to the positron flux (a 
few kiloparsecs). If the dark matter were locally dis- 
tributed completely evenly, the boost factor would be 
equal to one. Small scale clustering of dark matter, how- 
ever, enhances this quantity to a value expected to be in 
the range of roughly 2 to 5. Values much larger than this 
require large local clumps of dark matter and arc unlikely 
to be present [30]. 



IV. POSSIBLE VARIATIONS 

The results shown in the previous section may be 
modified with the choice of propagation parameters and 
Galactic dark matter distribution. In this section, we 
describe the effect of modifying these choices on the 
positron spectrum. 

The diffusion parameters which must be set in our for- 
malism are the diffusion zone width, 2L, the energy loss 
rate, b, and the diffusion constant, K, including its en- 
ergy dependence. The impact of varying these parame- 
ters is shown in figures 10-13. 

The effect of modifying the diffusion zone width is 
shown in figure 10. The flux at high energies is less af- 
fected by the choice of L than at lower energies. This 
is because distant annihilations produce positrons which 
are substantially degraded in energy before they reach 
Earth. Thus nearby annihilations produce most of the 
positrons at the high end of the spectrum. As the width 
of the diffusion zone is reduced, fewer positrons reach 
Earth at lower energies. 

Variations in the magnitude of the diffusion constant 
of Eq. 2 can also have an impact on the positron flux 
observed. As this quantity is increased, the lowest en- 
ergy positrons are most dramatically effected. Positrons 
of higher energies are actually mildly enhanced as the 
diffusion constant is increased. 
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FIG. 2. The spectrum of positrons, including the effects of 
propagation, from dark matter annihilations to b quark pairs. 
WIMP masses of 50, 100, 300 and 600 GeV were considered. 
A dark matter distribution with BF = 5 (see section III), 
p(local) = 0.43GeV/cm 3 and an annihilation cross section of 
av = 10 -25 cm 3 /s was used. The effects of solar modulation 
are not included. 
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FIG. 4. The spectrum of positrons, including the effects of 
propagation, from dark matter annihilations to e + e~ pairs. 
WIMP masses of 50, 100, 300 and 600 GeV were considered. 
A dark matter distribution with BF = 5 (see section III), 
p(local) = 0.43 GeV/cm 3 and an annihilation cross section of 
av = 10~ 26 cm 3 /s was used. The effects of solar modulation 
are not included. 
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FIG. 3. The spectrum of positrons, including the effects 
of propagation, from dark matter annihilations to ZZ and 
W + W~ pairs. WIMP masses of 100, 300 and 600 GeV were 
considered. Note the enhancement near mx/2 as compared to 
the case of annihilations to b quarks. A dark matter distribu- 
tion with BF = 5 (see section III), p(local) = 0.43GeV/cm 3 
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FIG. 5. The spectrum of positrons, including the effects of 
propagation, from dark matter annihilations to t + t~ pairs. 
WIMP masses of 50, 100, 300 and 600 GeV were considered. 
A dark matter distribution with BF = 5 (see section III), 
p(local) = 0.43 GeV/cm 3 and an annihilation cross section of 
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1 cm /s was used. The effects of solar modulation 



and an annihilation cross section of av = 10 J cm /s was are not included, 
used. The effects of solar modulation are not included. 
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FIG. 6. The positron fraction in the cosmic ray spectrum 
from dark matter annihilations to b quark pairs. WIMP 
masses of 50, 100, 300 and 600 GeV were considered. A 
dark matter distribution with BF — 5 (see section III), 
p(local) = 0.43GeV/cm 3 and an annihilation cross section 
of av = 10 -25 cm 3 /s was used. 
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FIG. 8. The positron fraction in the cosmic ray spec- 
trum from dark matter annihilations to e + e _ pairs. WIMP 
masses of 50, 100, 300 and 600 GeV were considered. A 
dark matter distribution with BF = 5 (see section III), 
p(local) = 0.43GeV/cm 3 and an annihilation cross section 
of av = 10 -26 cm 3 /s was used. 
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FIG. 7. The positron fraction in the cosmic ray spectrum 
from dark matter annihilations to ZZ and W + W~ pairs. 
WIMP masses of 100, 300 and 600 GeV were considered. 
A dark matter distribution with BF — 5 (see section III), 
p(local) = 0.43GeV/cm 3 and an annihilation cross section of 
av = 10~ 25 cm 3 /s was used. 
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FIG. 9. The positron fraction in the cosmic ray spec- 
trum from dark matter annihilations to t + t~ pairs. WIMP 
masses of 50, 100, 300 and 600 GeV were considered. A 
dark matter distribution with BF — 5 (see section III), 
p(local) = 0.43GeV/cm 3 and an annihilation cross section 
of av = 10~ 25 cm 3 /s was used. 
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Modifications of the positron energy loss rate (or equiv- 
alently the density of diffuse starlight and cosmic mi- 
crowave background photons) is tied to changes in the 
diffusion constant. Inspection of Eq. 1 reveals that the 
resulting spectral shape (up to normalization) is invariant 
with respect to the quantity K(E e +,x)/b(E e +,x). The 
normalization of the positron flux, on the other hand, 
scales inversely with either the diffusion constant or the 
energy loss rate. 

Finally, the choice of dark matter halo profile model 
considered can affect the positron flux observed. Thus 
far, our results have used an isothermal sphere dark mat- 
ter profile, described by 



p{r) = 0.43GeV/cm 



3 (2.8 kpc) 2 



.5kpc) 2 



(2.8 kpc) 2 



(5) 



where r is the distance from the Galactic Center. This 
profile exhibits a fairly flat behavior in the inner few kilo- 
parsecs of the Galaxy. Alternatively, we may consider a 
dark matter distribution favored by N-body simulations, 
such as the Navarro-Frenk- White (NFW) [31] profile: 



p(r) = 43 GeV/cm 3 (8-5/20) [1 + (8.5/20)] 2 
P{ ) u - 4J ^ cV / cm ( r /20kpc)[l + (r/20kpc)] 2 

or the Moore et al. profile [32]: 



. (6) 



p(r) = 0.43 GeV/cm 



: , (8.5/28) 1 - 5 [1 + (8.5/28) 1 - 5 ] 
(r/28 kpc) 1 - 5 [1 + (r/28 kpc) 1 - 5 ] ' 

(7) 



Unlike the isothermal sphere profile, the NFW and Moore 
et al. profiles have density cusps at the Galactic Center. 
Each of these profiles have been normalized to a local 
dark matter density of 0.43GeV/cm 3 . The only loca- 
tions which can deviate dramatically from this value are 
in the region near the Galactic Center. Sharply cusped 
profiles, the Moore et al. profile in particular, produce 
more positrons in the inner Galaxy and therefore more 
lower energy positrons are observed locally. This effect 
is shown in figure 13. The NFW and isothermal pro- 
files produce very similar results, with fewer low energy 
positrons than the Moore et al. case. Shown for com- 
parison is the result for a completely flat dark matter 
distribution, p(r) = constant = 0.43GeV/cm 3 . 



V. NEUTRALINO DARK MATTER 

In many models of supersymmetry, the Lightest Super- 
symmetric Particle (LSP) is a neutralino [33], a mixture 
of the superpartners of the photon, Z and neutral Higgs 
bosons. In models of supersymmetry with conserved R- 
parity, the LSP is stable and a potentially viable dark 
matter candidate. In this section, we will discuss the sig- 
natures of neutralino annihilations in the cosmic positron 
spectrum. 



The dominant annihilation channels of a neutralino 
depend on its composition. If the neutralino is mostly 
bino-like (the superpartner of the hypercharge gauge bo- 
son), it will annihilate largely to heavy fermions, i.e. 
bb, t + t~ and, if kincmatically allowed, ti. Annihilations 
to lighter fermions are rare as the low velocity neutralino 
annihilation cross section is chirality suppressed by a fac- 
tor of rn 2 /m 2 . Furthermore, the annihilation cross sec- 
tion to down-type fermions, such as bb, is enhanced in 
some diagrams by a factor of tan 2 f3. In much of the su- 
persymmetric parameter space, the LSP annihilates al- 
most entirely to bb with a small admixture of t~t+. If 
tan (3 is fairly small and the LSP mass is greater than the 
top mass, annihilations to ti can also be important. 

The lightest neutralino might not be mostly bino, how- 
ever. Alternatively, it may be largely higgsino-like (su- 
perpartner of the neutral Higgs bosons). In this case, 
the LSP may annihilate primarily to gauge boson pairs. 
The low velocity annihilation cross sections for a pure 
higgsino are approximately given by: 



< (TV >w+w- 



T F m W m x° 



* 1 7r(2m 2 - m 2 ,) 2 



and 



< (TV >z»Z°- 



m z 
m 2 „ 



G 2 F m%m^ Q 
27r(2m 2 - m|) 2 ' 



(8) 



(9) 



These cross sections can be quite large, reaching 3 x 
10~ 25 cm 3 /s at their maximum near m x o = 110 GeV. 
A mixed higgsino-bino LSP might also annihilate mostly 
to gauge bosons, but with a smaller cross section. 

As a third possibility, the lightest neutralino may also 
be a wino (the superpartner of the isospin gauge boson). 
Although wino-likc LSPs are uncommon in most super- 
symmetry breaking scenarios, in the interesting case of 
Anomaly Mediated Supersymmetry Breaking (AMSB), 
the lightest neutralino is a nearly pure wino. This is 
because the ratios of the gaugino masses do not follow 
the common hierarchy (Mi < < M3), but rather 
are proportional to /3-functions resulting in the ratios 
Mi : M 2 : M 3 = 2.8 : 1 : 7.1. When this ratio of 
Mi to M2 is inserted into the neutralino mass matrix, 
some generic phenomenological features of this model 
emerge. In addition to the LSP being a neutral wino, 
a chargino only a few hundred MeV heavier than the 
LSP is present. This can lead to efficient coannihilations 
in the relic density calculation. Much like the higgsino 
case, wino-like neutralinos in AMSB annihilate primarily 
via gauge bosons channels [34]. 
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FIG. 10. The effect of varying the diffusion zone width (2L) 
on the positron spectrum. Results for L = 2, 4 and 10 kpc 
are shown. In each case, WIMP annihilations are to b quark 
pairs and the WIMP's mass is 300 GeV. Other parameters 
are the same as in the previous figures. 
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FIG. 11. The effect of varying the diffusion constant on the 
positron spectrum. From top to bottom on the left side of the 
figure, results are shown for a diffusion constant 0.25 (dotted), 
0.5 (dashed), 1.0 (solid), 2.0 (dashed) and 4.0 (dotted) times 
the value shown in Eq. 2. In each case, WIMP annihilations 
are to b quark pairs and the WIMP's mass is 300 GeV. Other 
parameters are the same as in the previous figures. 
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FIG. 12. The effect of varying the energy loss rate on the 
positron spectrum. From top to bottom on the right side 
of the figure, results are shown for a rate 0.25 (dotted), 0.5 
(dashed), 1.0 (solid), 2.0 (dashed) and 4.0 (dotted) times the 
value shown in Eq. 3. In each case, WIMP annihilations are 
to b quark pairs and the WIMP's mass is 300 GeV. Other 
parameters are the same as in the previous figures. 
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FIG. 13. The effect of varying the dark matter halo profile 
on the positron spectrum. From top to bottom on the left 
side of the figure, results are shown for a Moore et al. pro- 
file (dotted), an Navarro- Frenk- White profile (dashed), and 
isothermal sphere profile (solid) and a completely flat distri- 
bution (dot-dashed). In each case, WIMP annihilations are 
to b quark pairs and the WIMP's mass is 300 GeV. Other 
parameters are the same as in the previous figures. 
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We now turn our attention to these various types of 
neutralino LSPs in the context of their annihilations pro- 
ducing positrons. The spectral shape of the positrons 
produced in the annihilation of neutralino LSPs is often 
dominated by either bb or gauge boson channels. In the 
bino case, as we said before, annihilations are typically 
almost entirely to bb with only a few percent admixture 
to t+t~. To see the effect of this admixture, we com- 
pare the positron spectrum from annihilations purely to 
bb with the spectrum from annihilations 90% to bb and 
10% to t+t~ . Below about E e + ~ m x o/10, the spectrum 
is not significantly effected by the tau component. At 
higher energies the tau channel becomes more important. 
For a 10% t + t" admixture, the flux at E e + ~ m x o /2 
is roughly doubled over the pure bb case. For an equal 
fraction of annihilations into bb and ti, the result is sim- 
ilar. These spectra are shown in figure 14. We have 
already shown the positron spectra for annihilations to 
gauge boson pairs, which are expected for higgsino-like 
or wino-likc ncutralinos. 
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FIG. 14. The spectrum of positrons, including the effects 
of propagation, from bino-like neutralino annihilations. The 
solid line is for annihilations purely to bb. The dashed line is 
for annihilations to bb with a 10% admixture to t + t~. The 
dotted line represents 50% to bb and 50% to ti. For each case 
a 300 GeV neutralino was considered. A dark matter distribu- 
tion with BF = 5 (see section III), p(local) = 0.43GeV/cm 
and an annihilation cross section of av — 10~ 25 cm 3 /s was 
used. The effects of solar modulation are not included. 

Predicting the spectral shape of the positrons produced 
in dark matter annihilations is only useful if the flux is 
large enough to be observed over background. The most 
important quantity in determining this flux is the dark 
matter annihilation rate, which in turn depends on the 
low velocity WIMP annihilation cross section. 

If the dark matter is made up thermal relics produced 
in the early universe, its density today can be calculated 
in a straight forward way. The result of this standard 



^ wimp h 2 = 



1.07 x 10 9 GeV~ 



x F 



fgZ{a + 3b/x F ) ' 



(10) 



where Mpi is the Plank mass, <?* is the number of rela- 
tivists degrees of freedom assessable at freeze-out and a 
and b are the first and second terms in the expansion of 
the annihilation cross section: < av >= a + bv 2 + 0(v 4 ). 
x F is the WIMP mass over the temperature at freeze-out, 
and is generally about 20 for a weakly interacting dark 
matter candidate. This expression neglects the possible 
role of coannihilations between the WIMP and other par- 
ticles during freeze-out. Coannihilations may be impor- 
tant when the LSP is only slightly lighter than the Next- 
to-Lightcst Supcrsymmetric Particle (NLSP). In such a 
case, the neutralino relic density can be depleted below 
the value of Eq. 10 and therefore a smaller annihilation 
cross section would be required to provide the measured 
relic abundance. 

For a weakly interacting dark matter candidate, Eq. 10 
numerically yields: 



WIMP 



2 3 x 10 7 cnrys 
(a + 36/20) 



(11) 



Matching this to the cold dark matter density observed 
by WMAP, tt W iMph 2 ~ 0.1, we arrive at 



a + 3&/20- 3 x 10" 26 cm a /s. 



(12) 



The low velocity annihilation cross section, a, for neu- 
tralino LSPs with a thermal relic density in the range 
of 0.095 < Q, x oh 2 < 0.129 (dark shading) and 0.06 < 
fl x oh 2 < 0.16 (light green shading) is shown in figures 15 
and 16 for models with tan/3 equal to 5 and 50, respec- 
tively. In the high tan (3 case shown in figure 16, low 
velocity annihilations to bb are highly efficient and, for 
neutralinos heavier than about 80 GeV, the low velocity 
cross section is often roughly the maximum value allowed 
by Eq. 12, < av >~ a~3x 10~ 26 cm 3 /s. For smaller 
values of tan (3, the tb-term in the velocity expansion of the 
annihilation cross section can be important, thus allow- 
ing a wide range of low velocity cross sections as shown 
in figure 15. 

The points shown in figures 15 and 16 were calculated 
using a variation of the DarkSusy package [24] , varying 
the parameters /i, M\, M2, M3, A t , and the sfermion 
masses randomly up to 10,000 GeV, and with either sign 
when appropriate. The pscudoscalar Higgs mass was var- 
ied up to 1 TeV. Although the complete MSSM contains 
more than 100 free parameters, the results shown in fig- 
ures 15 and 16 represent an illustrative sample of the 
possible models. Note that the sampling technique used 
in our monte carlo is somewhat limited for finding mod- 
els with an LSP heavier than several hundred GeV. Such 
models certainly exist, although they are probably also 
quite fine tuned. 
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FIG. 15. The neutralino annihilation cross section in the 
low velocity limit for parameters within the MSSM which 
yield a thermal relic density in the range measured by WMAP 
(^cdm/i 2 — 0.129 — 0.095, dark shading) or a somewhat larger 
range (JIcdm/i 2 = 0.16 — 0.06, light green shading). For all 
models shown, the ratio of Higgs expectation values, tan/3, is 
set to 5. 



o 10~ 





, , , , ' ' 


, , , , i , , 

j|0"3|"* ** 




26 3 


H§j #* * * * f M # J 1 
ElTfc^ ^ ****** 

*2r * * 

#§§ „ ***** 




"I** ***** 




27 _• 


wr * * *** 
B** ***** 

11 * 

St* * 
** 

* * 


** ** * • 
# * ** ** * 
* * 


* 


28 _ 


f. 

* 

» 






29 


. . . . i .... i 


. . . . i . . 





200 400 600 800 



m x o (GeV) 

FIG. 16. The neutralino annihilation cross section in the 
low velocity limit for parameters within the MSSM which 
yield a thermal relic density in the range measured by WMAP 
(Ocdm/i 2 — 0.129 — 0.095, dark shading) or a somewhat larger 
range (JIcdm/i 2 = 0.16 — 0.06, light green shading). For all 
models shown, the ratio of Higgs expectation values, tan/3, is 
set to 50. 

From the perspective of the calculation of the ther- 
mal relic density, higgsino or wino-like LSPs may appear 
disfavored, as their annihilation cross section is often too 
large to satisfy conditions such as those in Eq. 12. In fact, 
to provide the thermal relic density observed by WMAP 
a pure higgsino must have a mass of about 1 TeV. A pure- 
wino would need to be even more massive to provide the 
observed quantity of dark matter thermally. 

Alternatively, neutralino LSPs may be produced by 



non-thermal mechanisms thereby evading the require- 
ments of Eq. 12. In some scenarios, objects with small 
annihilation cross sections may freeze-out only to decay 
into LSPs at a later time. Long lived, but unstable states, 
such as gravitinos or Q-balls, would be examples of such 
intermediate states. In a non-thermal mechanism such 
as this, the LSP's annihilation cross section can be con- 
siderably larger than in the standard thermal relic: case 
[35]. 

For a more detailed discussion of neutralino dark 
matter in various supersymmtry breaking scenarios, see 
Ref. [36]. 

VI. KALUZA-KLEIN DARK MATTER 

In models with extra spatial dimensions, Kaluza-Klein 
(KK) excitations of Standard Model particles may ap- 
pear, and in some models may provide a viable dark 
matter candidate [37,38]. In particular, in models with 
Universal Extra Dimensions (UED), in which all of the 
fields of the Standard Model are free to propagate in 
the bulk, the Lightest Kaluza-Klein Particle (LKP) may 
be stable [39]. This is because the conservation of mo- 
mentum in higher dimensional space is conserved, and 
therefore the KK number of a state is conserved. For a 
UED model to be phenomenologically viable, the extra 
compact dimensions must be modded out by an orbifold. 
Orbifolding leads to terms which violate KK number, but 
may leave a conserved quantity, called KK-parity, which 
ensures that the LKP is stable in much the same that 
R-parity stabilizes the LSP in supersymmetry. 

At tree level, all of the first level KK modes have a 
mass corresponding to the compactifaction scale of the 
extra dimensions and their zero-mode mass: (m^l)i) 2 = 
(m t -0)i) 2 + 1/R 2 . Radiative corrections to the masses of 
the tree level spectrum break the (near) degeneracy of 
the first level KK states [40]. The lightest KK state is 
most naturally the level one excitation of the hypercharge 
gauge boson, B^l). For the remainder of this article, we 
will refer to this state simply as the LKP or Kaluza-Klein 
Dark Matter (KKDM). 

Although this LKP is similar to a bino-like neutralino 
in many phenomenological respects, its characteristics 
differ in significant ways. Perhaps most importantly for 
our purposes, the LKP is a boson and therefore its anni- 
hilations to fermions are not chirality suppressed. Thus, 
unlike neutralino dark matter, KKDM can annihilate di- 
rectly to e + e~, and t + t~, which each yield a gen- 
erous number of high energy positrons. The annihilation 
cross section of KKDM is nearly proportional to the hy- 
percharge of the final state fermions to the fourth power, 
therefore its annihilations are primarily to charged lep- 
tons (approximately 20% per generation) . The remaining 
annihilations are to up-type quarks (approximately 11% 
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per generation), neutrinos (approximately 1.2% per gen- 
eration), Higgs bosons (approximately 2.3%) and down- 
type quarks (approximately 0.7% per generation). 
The total annihilation cross section of KKDM is 

95g| 1.7 x 10" 26 cm 3 /s 

< aV >- 32W LKP " ml KP (TeV) ' {U > 

This annihilation cross section consists entirely of an a- 
term in the expansion < av >= a + bv 2 + 0(t> 4 ), i.e. 
b = and < av >= a. Comparing this expression with 
Eq. 12, we find that the thermal relic density of KKDM 
matches the value observed by WMAP for a mass of 
about rriLKP ~ 800 GeV, neglecting the effect of coan- 
nihilations [38]. Coannihilations can play an important 
role in the freeze-out of KKDM, however. It is somewhat 
expected that the loop-level radiative corrections to the 
KK spectrum should be much smaller than the tree-level 
contribution, thus there should be many KK states only 
slightly more massive than the LKP. Take, for example, 
one of the cases studied by Tait and Servant [38] . In this 
scenario, the first level KK excitations of ep_ were nearly 
degenerate with the LKP. The — annihilation 
cross section is only slightly larger than the LKP-LKP 
annihilation cross section and the e^J -LKP coannihila- 
tion cross section is substantially smaller. Thus each 
species freeze-out quasi-independently and after freeze- 
out, the remaining e^'s decay to LKPs thus enhancing 
the LKP relic density. This is in contrast to the coanni- 
hilation mechanism often considered for neutralinos, in 
which coannihilations between the LSP and other super- 
symmetric particles can be large and efficiently deplete 
the neutralino relic density. 

Tait and Servant find that for three generations of 
e^'s, each 1% heavier that the LKP, a mass oimLKP — 
550 — 650 GeV is needed to produce the observed density 
of LKP dark matter [38]. Furthermore, other coannihi- 
lation channels may certainly play an important role in 
determining the relic density of KKDM, plausibly low- 
ering the mass of the LKP considerably. Regardless of 
the coannihilation channels available, the LKP mass is 
constrained by electroweak precision measurements to be 
heavier than about 300 GeV [41]. 
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FIG. 17. The spectrum of positrons including the effects 
of propagation from Kaluza-Klein Dark Matter (KKDM) an- 
nihilations. Annihilations of KKDM produce equal fractions 
of t + t~ , fi + fj,~ and e + e~ pairs (approximately 20% each) as 
well as up-type quarks (approximately 11% per generation), 
neutrinos (approximately 1.2% per generation), Higgs bosons 
(approximately 2.3%) and down-type quarks (approximately 
0.7% per generation). Results for KKDM masses of 300 and 
600 GeV are shown. A dark matter distribution with BF = 5 
(see section III), /o(local) = 0.43GeV/cm 3 and an annihila- 
tion cross section of av = 10 -25 cm 3 /s was used. The effects 
of solar modulation are not included. 
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FIG. 18. The positron fraction in the cosmic ray spec- 
trum from Kaluza-Klein Dark Matter (KKDM) annihilations. 
KKDM masses of 300 and 600 GeV were considered. A 
dark matter distribution with BF — 5 (see section III), 
p(local) = 0.43GeV/cm and an annihilation cross section 
of av = 10 -25 cm 3 /s was used. 
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VII. THE HEAT RESULTS 

The HEAT (High-Energy Antimatter Telescopes) ex- 
periment's balloon flights in 1994-95 and 2000 have mea- 
sured the cosmic positron spectrum between energies of 
approximately 1 to 30 GeV [18-20]. When this data is 
presented as a positron fraction, an interesting feature 
appears. The origin of this feature, a bump appearing 
at 7-10 GeV, is not understood. Although many efforts 
have been made to fit this feature with a contribution 
from annihilating dark matter, it has been shown to be 
very difficult to satisfactorily match the spectrum of the 
bump [15-17,42]. 

Far more statistically significant than the bump- 
feature present in the HEAT data, however, is the over- 
all excess of positrons above the background prediction 
at energies above about 7 GeV. This, in principle, can 
certainly be accommodated by a contribution to the 
positron spectrum from dark matter annihilations, al- 
though rather large cross sections or boost factors are 
required to do so. 

To judge the degree that adding a dark matter com- 
ponent to the positron spectrum improves the fit to the 
HEAT data, we calculate the x 2 , 

^2 _ (^Obs — VpM — -^BG) 2 

AjVobs 

where the sum is over energy bins (one for each HEAT 
error bar), Aobs is the number of events observed in that 
bin, Adm is the number of events predicted from the 
annihilating dark matter contribution, Nbg is the num- 
ber of events predicted from the background contribution 
and A jv 0b3 is the error associated with the measurement 
in each energy bin. 

The \ 2 we fi n d f° r t ne background-only curve is 47.2 
over 12 degrees of freedom (error bars). This x 2 of al- 
most 4 per degree of freedom represents a very poor fit 
to the HEAT data. The overall quality of this fit can 
be dramatically improved by including a new component 
from dark matter annihilations. 

Considering first annihilations primarily to bb, we find 
that only rather heavy WIMPs (to > 200 GeV) can ac- 
commodate the HEAT data. This can be seen from 
the 100 GeV WIMP curve shown in figure 19. For this 
curve, the flux at low energies consistently exceeds the 
HEAT measurements while it falls below the measure- 
ments above 10 GeV. We find that for the best-fit nor- 
malization (av = 9x 10 _26 cm 3 /s, BF = 5, pi O cai=0.43 
GeV/cm 3 ) a x 2 of 23.5, or almost 2 per degree of free- 
dom, is obtained, representing a rather poor fit. For a 
300 or 600 GeV WIMP annihilating to bb, however, we 
find a x 2 of 13.6 and 9.7, respectively. In these cases, the 
normalization must be increased significantly, however. 
Using BF = 5 and pi O cai=0.43 GeV/cm 3 , we must re- 
quire av = 2.5 x 10~ 25 cm 3 /s and av = 5.5 x 10 -25 cm 3 /s 



to achieve these good fits. These cross sections are each 
well above the value acceptable for a thermal relic (see 
Eq. 12). 
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FIG. 19. The positron fraction for dark matter annihila- 
tions to bb for WIMP masses of 100, 300 and 600 GeV (bot- 
tom to top on the right side of the figure). In each case, the 
normalization of the dark matter contribution was chosen to 
maximize the quality of the fit to the HEAT data (shown as 
error bars). We find x 2 ' s P er degree of freedom of 23.5/12, 
13.6/12 and 9.7/12 for 100, 300 and 600 GeV WIMPs, respec- 
tively. Considering a dark matter distribution with BF — 5 
(see section III) and p(local)=0.43 GeV/cm 3 , annihilation 
cross sections of av = 9 x 10" 26 cm 3 /s, 2.5 x 10~ 25 cm 3 /s 
and 5.5 x 10~ 25 cm 3 /s are needed for these fits for 100, 300 
and 600 GeV WIMPs, respectively. The solid line is the back- 
ground-only prediction. 

It is worth noting that the situation shown in figure 19 
can be improved if more favorable diffusion and halo pa- 
rameters are used. For example, if we use a diffusion 
constant twice the size of our standard choice (shown in 
Eq. 2), we can improve the fit to the HEAT data sub- 
stantially, consistent with WIMPs as light as 30-40 GeV. 
These results are shown in figure 20. 

Next, considering annihilations primarily to gauge 
bosons, such as would be predicted for a higgsino or wino- 
like neutralino, we find that a reasonable fit to the HEAT 
data can be achieved for a WIMP with any mass above 
the gauge boson pair production threshold (although the 
fit improves for a heavier WIMP). In this case, we do 
not require the large diffusion constant used in figure 20. 
These results are shown in figure 21. 
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FIG. 20. The same as in figure 19, except with a diffusion 
constant twice as large. WIMPs ol mass 50, 100 and 300 GeV 
are shown (100, 300, 50 from top to bottom on the right). 
Here we find x 2 ' s P er degree of freedom of 7.5/12, 9.7/12 and 
9.5/12 for 50, 100 and 300 GeV masses, respectively. Consid- 
ering a dark matter distribution with BF = 5 (see section III) 
and p(local)=0.43 GeV/cm 3 , annihilation cross sections of 
av = 6 x 10~ 26 cm 3 /s, 1.5 x 10 -26 cm 3 /s and 4.5 x 10 -25 cm 3 /s 
are needed for these fits for 50, 100 and 300 GeV WIMPs, re- 
spectively. The solid line is the background-only prediction. 
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FIG. 21. The positron fraction for dark matter annihila- 
tions to ZZ and W + W~ pairs for WIMP masses of 300, 600 
and 100 GeV (bottom to top at 30 GeV). In each case, the nor- 
malization of the dark matter contribution was chosen to max- 
imize the quality of the fit to the HEAT data (shown as error 
bars). We find x 2 ' s P er degree of freedom of 15.3/12, 11.4/12 
and 8.8/12 for 100, 300 and 600 GeV WIMPs, respectively. 
Considering a dark matter distribution with BF = 5 (see 
section III) and p(local)=0.43 GeV/cm 3 , annihilation cross 
sections of av = 1.4 x 10~ 25 cm 3 /s, 3.6 x 10~ 25 cm 3 /s and 
7.6 x 10~ 25 cm 3 /s are needed for these fits for 100, 300 and 
600 GeV WIMPs, respectively. The solid line is the back- 
ground-only prediction. 
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FIG. 22. The positron fraction for Kaluza-Klein Dark Mat- 
ter (KKDM) annihilations for WIMP masses of 300 and 600 
GeV (top to bottom on the right side of the figure). In each 
case, the normalization of the dark matter contribution was 
chosen to maximize the quality of the fit to the HEAT data 
(shown as error bars). We find x 2 's per degree of freedom of 
10.4/12 and 9.2/12 for 300 and 600 GeV WIMPs, respectively. 
Considering a dark matter distribution with BF — 5 (see sec- 
tion III) and p(local)=0.43 GeV/cm 3 , annihilation cross sec- 
tions of av = 5x 10~ 25 cm 3 /s and 1.2 xlO -24 cm 3 /s are needed 
for these fits for 300 and 600 GeV WIMPs, respectively. The 
solid line is the background-only prediction. 

As with the case of annihilations to bb, the cross section 
needed for annihilations to gauge bosons to fit the HEAT 
data are quite large, 1-8 x 10 -25 cm 3 /s for the mass range 
we studied with a boost factor of 5 and a local dark 
matter density of 0.43 GeV/cm 3 . Although a thermal 
relic cannot have such a large annihilation cross section, 
it is interesting to note that this is roughly the range 
predicted for a neutralino LSP in the Anomaly Mediated 
Supersymmetry Breaking (AMSB) scenario with a non- 
thermal production mechanism. 

Lastly, we show the fit of the positron spectrum from 
KKDM annihilations to the HEAT results in figure 22. 
For both masses we consider (300 and 600 GeV) , we find 
an excellent fit to the data [17]. Again, the annihilation 
cross sections needed are quite large. 

VIII. SENSITIVITY OF PAMELA AND AMS-02 

Given the difficulty of a thermal dark matter relic to 
produce the positron flux observed by HEAT and the 
uncertainties regarding the effects of solar modulation at 
GeV energies, it is certainly possible that HEAT has not 
revealed a signature of dark matter. To pursue this ques- 
tion further, future experiments will be required to mea- 
sure the positron spectrum with greater precision and at 
higher energies. The PAMELA and AMS-02 experiments 
will have the ability to perform such measurements. 

In this section, we will not assume a positron spectrum 
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consistent with the measurements of HEAT, but rather 
attempt to determine at what dark matter annihilation 
rate (or cross section) such a signature could be identified 
in future experiments. 

At a minimum, for an experiment to claim the ob- 
servation of a signature of dark matter annihilation, the 
spectral data must be statistically distinct from the pre- 
dicted background spectrum. To evaluate this, we can 
calculate the x 2 of a given model for a set of data and 
compare it to the background only case. 

To determine the x 2 01 a se t of data over the expected 
background, we perform the following sum: 



X 



E 



(N ohs - N BG ) 2 
Nobs 



(15) 



where the sum is over energy bins, Nq^ s is the number 
of events observed in that bin and -/Vbg is the number 
of events predicted from the background contribution. 
Here we have assumed Gaussian errors. We have chosen 
to adopt energy bins of width A(logE') = 0.60 below 40 
GeV and A(log£) = 0.66 above 40 GeV. 

Just as we can convert the fluxes of figures 2-5 to 
positron fractions by using fluxes of background positrons 
and electrons, we can replace the fluxes (or numbers 
of events) used in Eq. 15 with ratios of positrons to 
positrons+electrons to reduce the effects of solar mod- 
ulation in our results. Because there are substantially 
more electrons than positrons observed, we can assume 
that there are negligible errors associated with the elec- 
tron flux. 

For concreteness, we first consider the specific example 
of a 100 GeV WIMP annihilating to gauge bosons with a 
cross section of 10 -25 cm 3 /s, a boost factor of 5 and a lo- 
cal dark matter density of 0.43 GeV/cm 3 . The positron 
fraction found for this scenario is shown in figure 23 along 
with the error bars after three years of observation by 
PAMELA (light-red) and AMS-02 (dark-blue). Follow- 
ing Eq. 14, the total \ 2 obtained for this flux is 1045.3 
for PAMELA and 22,944 for AMS=02, both clearly dis- 
tinguishable from the background. 

If the annihilation rate is smaller, identifying such a 
feature will be more difficult. If we consider an annihi- 
lation rate ten times smaller than that used in figure 23, 
we find the results shown in figure 24. It is clear to the 
eye that the excess over background in this case is far 
less statistically significant. With three years of observa- 
tion by PAMELA, the x 2 obtained is only 15.1. AMS-02, 
however, still obtains a x 2 of 1045.3. 
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FIG. 23. The positron fraction for dark matter annihila- 
tions to gauge bosons with an annihilation cross section of 
10~ 25 cm 3 /s, a boost factor (see section III) of 5, a local dark 
matter density of 0.43 GeV/cm 3 and a WIMP mass of 100 
GeV. The light-red and dark-blue error bars shown are those 
projected for the PAMELA and AMS-02 experiments, respec- 
tively; each after three years of observations. The solid line 
is the background-only prediction. 
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FIG. 24. The positron fraction for dark matter annihila- 
tions with a WIMP with the same characteristics as in fig- 
ure 23, except with an annihilation cross section ten times 
smaller, 10 -26 cm 3 /s. Again, the light-red and dark-blue error 
bars shown are those projected for the PAMELA and AMS-02 
experiments, respectively; each after three years of observa- 
tions. The solid line is the background-only prediction. 
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FIG. 25. The positron fraction for dark matter annihila- 
tions to bb with an annihilation cross section of 10~ 25 cm 3 /s, 
a boost factor (see section III) of 5, a local dark matter den- 
sity of 0.43 GeV/cm 3 and a WIMP mass of 100 GeV. The 
light-red and dark-blue error bars shown are those projected 
for the PAMELA and AMS-02 experiments, respectively; each 
after three years of observations. The solid line is the back- 
ground-only prediction. 
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FIG. 27. The positron fraction for dark matter anni- 
hilations to e + e~ with an annihilation cross section of 
10~ 26 cm 3 /s, a boost factor (see section III) of 5, a local dark 
matter density of 0.43 GeV/cm 3 and a WIMP mass of 100 
GeV. The light-red and dark-blue error bars shown are those 
projected for the PAMELA and AMS-02 experiments, respec- 
tively; each after three years of observations. The solid line 
is the background-only prediction. 
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FIG. 26. The positron fraction for dark matter annihila- 
tions with a WIMP with the same characteristics as in fig- 
ure 25, except with an annihilation cross section ten times 
smaller, 10 -26 cm 3 /s. Again, the light-red and dark-blue error 
bars shown are those projected for the PAMELA and AMS-02 
experiments, respectively; each after three years of observa- 
tions. The solid line is the background-only prediction. 
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FIG. 28. The positron fraction for dark matter annihila- 
tions with a WIMP with the same characteristics as in fig- 
ure 27, except with an annihilation cross section ten times 
smaller, 10 -27 cm 3 /s. Again, the light-red and dark-blue error 
bars shown are those projected for the PAMELA and AMS-02 
experiments, respectively; each after three years of observa- 
tions. The solid line is the background-only prediction. 
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The situation is somewhat less optimistic for dark mat- 
ter particles which annihilate to b quarks. These results 
are shown in figures 25 and 26. Again, with the large 
cross section of 10 -25 cm 3 /s (assuming BF = 5, p = 0.43 
GeV/cm 3 and a 100 GeV WIMP), there is a significant 
excess over the expected background (see figure 25). We 
calculate x 2 's of 9439 and 430 for AMS-02 and PAMELA, 
respectively. For the smaller cross section of 10~ 26 cm 3 /s, 
the excess becomes very small and the \ 2 falls to 5.4 and 
118 at PAMELA and AMS-02 (see figure 26), impossible 
to resolve with PAMELA and a challenge for AMS-02. 

The x 2 °f the deviation from the background predic- 
tion is not the only measure of an experiment's ability to 
identify a contribution from dark matter annihilations. 
This becomes apparent when considering the case of dark 
matter which annihilates directly to positrons, i.e. e + e~ 
pairs. The results for this case is shown in figures 27 
for a cross section of 10 _26 cm 3 /s and in figure 28 with 
10~ 27 cm 3 /s. Corresponding to figure 27, we find x 2 ' s 
of 17,529 and 799 for AMS-02 and PAMELA, respec- 
tively. For the smaller cross section of figure 28, we find 
X 2 's of 334 and 15.2. Note that these x 2 ' s are some- 
what smaller than those corresponding to annihilations 
to gauge bosons (figures 23 and 24). Visually comparing 
these cases, it is obvious that the sudden discontinuity 
present in the e + e~ case, but absent in the spectrum 
from annihilations to gauge bosons (or bb, etc.), provides 
a clear signature of a new positron contribution which 
could not be incorporated into a modified background 
model. Despite the somewhat smaller \ 2 values of the 
e + e~ case, this spectrum has distinctive features which 
more than make up for this difference. 

Lastly, our results for Kaluza-Klein Dark Matter 
(KKDM) are shown in figures 29 and 30. We calculate 
X 2 values of 393 and 8626 for the annihilation rate shown 
in figure 29, for the PAMELA and AMS-02 experiments, 
respectively. For the smaller annihilation rate used in 
figure 30, we find x 2 's of 5.9 and 129 for the two experi- 
ments. Note that in these two figures a mass of 300 GeV, 
corresponding to the minimum KKDM mass allowed by 
electroweak precision measurements, is used rather than 
the 100 GeV mass used in figures 23-28. 
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FIG. 29. The positron fraction for Kaluza-Klein Dark Mat- 
ter (KKDM) annihilations. A WIMP mass of 300 GeV, the 
minimum mass allowed for KKDM by electroweak precision 
measurements, was used. The annihilation rate shown corre- 
sponds to a boost factor (see section III) of 5, a local dark 
matter density of 0.43 GeV/cm 3 and an annihilation cross sec- 
tion of 10 _25 cm 3 /s. The light-red and dark-blue error bars 
shown are those projected for the PAMELA and AMS-02 ex- 
periments, respectively; each after three years of observations. 
The solid line is the background-only prediction. 
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FIG. 30. The positron fraction for dark matter annihila- 
tions with a WIMP with the same characteristics as in fig- 
ure 29, except with an annihilation rate ten times smaller. 
Again, the light-red and dark-blue error bars shown are those 
projected for the PAMELA and AMS-02 experiments, respec- 
tively; each after three years of observations. The solid line 
is the background-only prediction. 
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IX. REACH OF PAMELA AND AMS-02 

To quantitatively assess the reach of PAMELA and 
AMS-02, we calculate the annihilation rate needed to be 
distinguished from the background at the 95% confidence 
level. The 95% confidence level corresponds to a x 2 per 
degree of freedom of approximately 1.3. For our method 
of energy binning, this corresponds to a total x 2 of about 
30. Although this value can vary somewhat depending 
on the energy bins used and other assumptions, our over 
results are only slightly modified by these variations. 

For a given WIMP mass and annihilation mode(s), we 
can find the dark matter annihilation rate that corre- 
sponds to the 95% confidence level sensitivity. The anni- 
hilation rate can be further separated into the low veloc- 
ity annihilation cross section and the Boost Factor (BF) 
(see section III). Throughout this section, we assume 
a mean local dark matter of p = 0.43 GeV/cm 3 . The 
annihilation rates which are observable in three years of 
observation by PAMELA or AMS-02 are shown in fig- 
ures 31, 32 and 33. 

For a bino-like, neutralino dark matter candidate, we 
see the range of annihilation cross sections which are as- 
sessable to PAMELA and AMS-02 in figure 31. With 
typical boost factors of a few, PAMELA may be capa- 
ble of observing a bino-like neutralino with a cross sec- 
tion near the maximum allowed for a thermal relic (a 
few times 10 _26 cm 3 /s). The prospects for detection are 
particularly good for models with large values of tan/3 
(see figures 15 and 16). For AMS-02, the sensitivity 
is considerably improved, reaching cross sections around 
10~ 27 cm 3 /s. We find very similar results if we consider 
a WIMP_which annihilates 90% to bb and 10% t+t" or 
50% to bb and 50% ti (see section V). 

For a higgsino or wino-like LSP, our results are shown 
in figure 32. Also shown in this figure are the predicted 
annihilation cross sections to gauge boson pairs for a pure 
higgsino and a pure wino within the Anomaly Mediated 
Supersymmetry Breaking (AMSB) scenario. Due to the 
nearly degenerate chargino present in the AMSB sce- 
nario, LSP annihilation cross sections to gauge bosons 
are very large, as can be seen in the figure. After three 
years, PAMELA will be able of testing the AMSB sce- 
nario with masses as large as 550 GeV and 1 TeV, for 
boost factors in the range of 1 to 5. AMS-02 will be sen- 
sitive to this scenario for masses up to 2-3 TeV. For a 
higgsino, PAMELA will be sensitive up to 230-380 GeV 
and while AMS-02 can probe as high as 400-650 GeV. If 
a mixed gaugino-higgsino is the LSP, results somewhere 
between those shown in figures 31 and 32 will apply. 



10 



-24 



10 



-28 



PAMELA, BF=1 




800 



200 400 600 

mwiMP (GeV) 

FIG. 31. The ability of PAMELA and AMS-02 after three 
years of observation to detect positrons from dark matter an- 
nihilations into bb at the 95% confidence level. Shown from 
top to bottom are the reach of PAMELA with a boost fac- 
tor of 1 (dashed), PAMELA with a boost factor of 5 (solid) 
and AMS-02 with a boost factor of 5 (solid). Not shown is 
the reach of AMS-02 to dark matter with a boost factor of 1, 
which contour falls nearly on top of the PAMELA case for a 
boost factor of 5. Note that the maximum low velocity an- 
nihilation cross section for a thermal relic is approximately 
3 x 10" 26 cm 3 /s. 
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FIG. 32. As in figure 31, except for dark matter which an- 
nihilates into gauge bosons. Also shown as thick solid lines 
are the predicted annihilation cross sections for a pure hig- 
gsino-LSP (lower, blue) and a nearly pure wino-LSP in the 
AMSB scenario (upper, green). See section V for more de- 
tails. 
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FIG. 33. As in figure 31, except for Kaluza-Klein Dark 
Matter (KKDM), annihilating in the modes described in sec- 
tion VI. Shown as a solid blue line is the annihilation cross 
section predicted for KKDM. PAMELA will be able to test for 
the presence of KKDM up to masses of 350-550 GeV, while 
AMS-02 will test up to masses of 550-1000 GeV. 

Our results for Kaluza-Klein Dark Matter (KKDM) 
are shown in figure 33. We find that PAMELA will test 
this scenario for masses up to 350-550. AMS-02 will reach 
up to 1 TeV. 

X. COMPARISONS TO OTHER DETECTION 
METHODS 

In this section, we will briefly compare the prospects 
for observing particle dark matter with cosmic positrons 
to dark matter detection using other techniques. 

The prospects for direct detection experiments [6], 
which look for the effects of WIMPs elastically scatter- 
ing off of an detector, are difficult to compare to those 
of cosmic positron experiments. This is because direct 
detection experiments largely depend on the elastic scat- 
tering cross section of a WIMP with the target material, 
whereas the success of positron experiments depends on 
the WIMP's annihilation cross section (at low velocities). 
These two cross sections are not directly determined by 
each other, and are very difficult to compare in any gen- 
eral way. 

The prospects for dark matter detection with neutrino 
telescopes (by observing neutrinos produced in the anni- 
hilations of WIMPs trapped in the Sun or Earth) [9] are 
also difficult to compare to positron experiments. Again, 
these experiments largely depend on the elastic scatter- 
ing cross section of a WIMP with nucleons rather than a 
WIMP's annihilation rate. 

Searches for gamma-rays produced in WIMP annihi- 
lations [8] is the method of indirect detection which has 
been studied in the greatest detail. The prospects for 
this method depend strongly on the distribution of dark 
matter in our Galaxy, particularly in the Galactic Cen- 



ter. If a high density of dark matter (a cusp or a spike) 
is present in the Galactic Center, it has been shown that 
future gamma-ray experiments such as GLAST or HESS 
would likely be able to detect the annihilation signal. 
Very recently, however, HESS has announced the detec- 
tion of a bright gamma-ray source at the Galactic Center 
which may interfere with the future detection of gamma- 
rays produced in dark matter annihilations [43]. This 
HESS source has a spectrum consistent with a power 
law extending at least to several TeV. Although it has 
been shown that this, in principle, could be the product 
of dark matter annihilation, it would require an exceed- 
ingly heavy WIMP (above 12 TeV), which is probably 
less likely than the presence of a astrophysical source at 
this location. With such a source coincident with the 
dark matter annihilation region at the Galactic Center, 
the continuum component of the dark matter gamma-ray 
flux will be difficult to identify, leaving only the much 
more faint line emission to be observed in future gamma- 
ray experiments. 

In addition to positrons, other anti-matter species can 
be used to search for evidence of dark matter annihila- 
tions. In particular, observations of cosmic anti-protons 
[10,11] and anti-deuterons [11,12] may be capable of re- 
vealing signatures of particle dark matter. Measure- 
ments of the cosmic anti-proton spectrum have excellent 
prospects for dark matter detection, although they de- 
pend critically on unknown parameters of the Galactic 
halo, so are difficult to assess with much confidence. In 
particular, the width and extent of the galactic diffu- 
sion zone can dramatically alter the rates and spectra 
predicted for anti-protons. Positrons, on the other hand, 
lose energy much more rapidly and therefore sample only 
the surrounding few kpc of the dark matter halo, thus 
making the diffusion zone boundary conditions less im- 
portant. Anti-deuteron studies certainly look promising, 
although less attention has been given to these than to 
other detection channels thus far. 

XI. CONCLUSIONS 

Future measurements of the cosmic positron spectrum 
represent one of the most promising methods for detect- 
ing particle dark matter. The PAMELA and AMS-02 ex- 
periments each represent major steps forward in the pre- 
cision measurement of the cosmic positron spectrum, for 
the first time measuring up to hundreds of GeV. If dark 
matter consists of weakly interacting particles present in 
our Galactic halo, it is fairly likely that these experiments 
will be sensitive to it. 

If the dark matter is a supersymmetric neutralino, the 
prospects for its detection by PAMELA and AMS-02 de- 
pend on its composition. If it is a nearly pure higgsino 
(superpartner of the neutral Higgs bosons), PAMELA 
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will be sensitive to masses up to 230 to 380 GeV, af- 
ter three years of observation, depending on the degree 
of local inhomogeneities in the dark matter distribution. 
AMS-02 will be sensitive up to 400 to 650 GeV for a hig- 
gsino LSP. If the LSP is bino-like, as is the case in much 
of the parameter space of the Minimal Supersymmetric 
Standard Model (MSSM), the prospects for these experi- 
ments depend strongly on the WIMP's annihilation cross 
section (at low velocities). PAMELA will be sensitive 
to cross sections up to the order of a few times 10~ 26 
cm 3 /s while AMS-02 can improve on this by a factor of 
about four. Supersymmetric models with large values 
of tan [3 are particularly likely to fall near or above this 
range. Thirdly, if supersymmetry comes in the form of 
Anomaly Mediated Supersymetry Breaking (AMSB), the 
wino LSP will be quite observable by both PAMELA and 
AMS-02, probing such models up to multi-TcV masses. 

If the dark matter consists of Kaluza-Klein excitations 
of Standard Model fields (the excitation of the hyper- 
charge gauge boson, in particular), their annihilations 
should produce an observable positron flux to AMS-02, 
and possibly to PAMELA. PAMELA will be sensitive to 
KKDM masses up to 350 to 550 GeV, depending on the 
degree of local inhomogeneities in the dark matter dis- 
tribution. AMS-02 will be sensitive up to 600 GeV to 1 
TeV. 

The prospects for the detection of WIMPs from their 
positron signatures with PAMELA and AMS-02 are quite 
encouraging. We look forward to the deployment of these 
experiments in the coming years. 
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